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Abstract 

The rate of dd^i muonic molecule resonant formation in d^ atom collision with a condensed 
deuterium target is expressed in terms of a single-particle response function. In particular, ddfi 
formation in solid deuterium at low pressures is considered. Numerical calculations of the rate in 
the case of fee polycrystalline deuterium at 3 K have been performed using the isotropic Debye 
model of solid. It is shown that the energy-dependent ddfi formation rates in the solid differ 
strongly from those obtained for D2 gaseous targets, even at high d/i kinetic energies. Monte 
Carlo neutron spectra from dd fusion in dd^i molecules have been obtained for solid targets with 
different concentrations of ortho- and para-deuterium. The recent experimental results performed 
in low pressure solid targets (statistical mixture of ortho-D2 and para-D2) are explained by the 
presence of strong recoilless resonance peaks in the vicinity of 2 meV and very slow deceleration 
of dn atoms below 10 meV. A good agreement between the calculated and experimental spectra 
is achieved when a broadening of D2 rovibrational levels in solid deuterium is taken into account. 
It has been shown that resonant dd^i formation with simultaneous phonon creation in solid gives 
only about 10% contribution to the fusion neutron yield. The neutron time spectra calculated for 
pure ortho-D2 and para-D2 targets are very similar. A practically constant value of the mean dd^ 
formation rate, observed for different experimental conditions, is ascribed to the fact that all the 
recent measurements have been performed at temperatures T < 19 K, much lower than the target 
Debye temperature 0d ~ 110 K. In result, the formation rate, obtained in the limit T/Bd ^ 1, 
depends weakly on the temperature. 
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I. INTRODUCTION 



Theoretical study of resonant formation of the muonic molecule ddfi in condensed deu- 
terium targets is the main subject of this paper. The resonant ddfi formation, first observed 
by Dzhelepov and co-workers 0], is a key process of muon catalyzed fusion (/iCF) in deu- 
terium (see e.g. reviews 0, A muonic deuterium atom dfi is created when a negative 
muon /x~ is captured into an atomic orbital in a deuterium target. After dn deexcitation 
to the 15* state and slowing down, the ddfi molecule can be formed in d/i atom collision 
with one of the D2 target molecules. The resonant formation is possible due to presence 
of a loosely bound state of dd^, characterized by the rotational number J = 1 and vibra- 
tional number v = 1, with binding energy ^ 1.97 eV. This energy, according to 
the Vesman mechanism [Q, is completely transferred to excited rovibrational states of the 
molecular complex [{ddfi)dee]. The scheme of calculation of dd^ formation rate in gaseous 

|], and has lead to a good agreement 
]lO| . On the other hand, this theory, 



deuterium has been developed for many years ^ 
with the experiments performed in gaseous targets 
when directly applied to solid deuterium targets, leads to strong disagreement with the ex- 
perimental results in], |l 
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Therefore, it is necessary to calculate the ddfi formation rate 
with solid state effects taken into account, which is the main purpose of this paper. 

Our calculations are based on the theoretical results (transition matrix elements, reso- 
nance energies) obtained in the case of ddfi formation in a single D2 molecule. In Sec. II 
the main formulas used for this case are briefly reported. A general formula for the energy- 
dependent ddfi formation rate in a D2 condensed target is derived in Sec. |T|, using the Van 
Hove formalism of the single-particle response function [|14|. This formula is then applied 



(Sec. [TVD for harmonic solid targets, in particular for a cubic Bravais lattice. A phonon 
expansion of the response function is used to study phonon contributions to the resonant 
formation. Numerical results for 3 K zero pressure frozen deuterium targets (TRIUMF ex- 
perimental conditions [|1T], |13|)) with the fee polycrystalline structure, are shown in Sec. |V]. 



The formation rates have been calculated assuming the isotropic Debye model of the solid 
and the values of Debye temperature and lattice constant observed in neutron scattering 
experiments. 

The calculated rates of resonant ddfi formation and back decay have been used for Monte 
Carlo simulations of dd fusion neutron and proton time spectra. Since the initial distributions 
of IS" muonic atom energy contain contributions from hot d/i's (~ 1 eV) |jl5, 
of slow deceleration of dfi atoms below 10 meV 
in Sec. 



161, influence 



FFI on these time spectra is investigated 



VI 



The simulations take into account the processes of incoherent and coherent dfi 
atom scattering in solid deuterium. In particular, the Bragg scattering, phonon scattering, 
and rovibrational transitions in D2 molecules are included. We consider a dependence of 
the resonant formation rate and time spectra on broadening of the rovibrational D2 energy 
levels, due to the binding of the molecules in the lattice HI 



Since it has been predicted in Refs. |jT9|, ^ that strong ddfi formation takes place only 



in solid para-D2, study of this process in pure ortho-D2 and para-D2 targets is another aim 
of this work. The neutron spectra calculated for these two solids are discussed in Sec. [V^. 
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II. RESONANT FORMATION IN A FREE MOLECULE 



First we consider resonant formation of the ddfi molecule in the following reaction 

(d/i)^ + [{ddf^Ys^dee] ^^^^ , (1) 

where D2 is a free deuterium molecule in the initial rovibrational state {viKi) and the 
total nuclear spin I. The muonic atom d^ has total spin F and CMS kinetic energy e. 
The complex [{dd^)dee\ is created in the rovibrational state [vfKf) and the molecular ion 
dd/i, which plays the role of a heavy nucleus of the complex, has total spin S. The rate 
^l^Ki,ufKf process above depends on the elastic width ^l^Kf,ViKi [{ddij)dee] complex 

decay |2^, in reactions 

{{ddii)i^dee\^^^^ (2) 

— s> stabilization processes, 

h 

where A/ is the total rate of the stabilization processes, i.e. deexcitation and nuclear fusion 
in ddii 

'/i + t + P + 4.0 MeV 
(idyu ^ + ^He + n + 3.3 MeV (3) 
/i3He + n + 3.3 MeV. 

When fusion takes place, the muon is generally released and can again begin the /iCF cycle. 
However, sometimes the muon is captured into an atomic orbital of helium (sticking), which 
stops further reactions. 

The value of ^ufxf UiKi given in atomic units (e = ^ = mg = 1) by the formula 

/d^k 
(4) 

where Vif{e) is a transition matrix element and £jj is a resonance energy defined in Ref. [Q. 
The factor Aij is due to averaging over initial and summing over final projections of spins 
and angular momenta of the system. Vector k is the momentum of relative d^ and D2 
motion 

e = e/2M , (5) 
and Ai is the reduced mass of the system. Integration of Eq. over k leads to 



TT 



Since ^ufKf,UiKi niuch lower (~ 10 ^ meV) than e, Vesman's model can be applied 

and the energy-dependent resonant formation rate has the Dirac delta function profile 



= 2nNB,j\V,fie)\^5ie - e,f) . (7) 
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where is the density of deuterium nuclei in the target. According to Ref. P] the coefficients 
Aif and Bif in the above equations are equal to 



where 



(8) 



Ws, = i2F + l)[iy^ 

m) = 



2 



I for Ki = 0, 
I for ir, = 1 , 



(9) 



and the curly brackets stand for the Wigner 6j symbol. In formula the usual Boltzmann 
factor describing the population of rotational states in a gas target is omitted because we 
calculate the formation rate separately for each initial rotational state. If the muonic atoms 
in a gas have a steady kinetic energy distribution f{e,T) at target temperature T, Eq. (0) 
can be averaged over the atom motion leading to a mean resonant rate A^^. jy^i^^(T). 



III. RESONANT FORMATION IN A CONDENSED TARGET 

Since a muonic deuterium atom can be approximately treated as a small neutron-like 
particle, methods used for description of neutron scattering and absorption in condensed 
matter are applicable in the case ddfi formation in dense deuterium targets. Below we adapt 
the method developed by Lamb |2^, and then generalized by Singwi and Sjolander 



using the Van Hove formalism of the single-particle response function Si [|T4| , for calculation 
of the resonant dd/j, formation rates. 

A Hamiltonian iftot of a system, consisting of a dfi atom in the IS* state and a heavy 
condensed D2 target, can be written down as follows 

^tot = ^V^^^ + H,,{r,) + Hj,M + V{n, Q„ Q,) + H, (10) 

where M^^ is the dfi mass and R^^ denotes the position of dfi center of mass in the coor- 
dinate frame connected with the target (see Fig. |l|). Operator Hdfj, is the Hamiltonian of 
a free dfi atom, ri is dn internal vector; H^^ denotes the internal Hamiltonian of a free D2 
molecule. It is assumed that ddn formation takes place in collision with the 1-th D2 target 
molecule. The position of its mass center in the target frame is denoted by R^; Qi is a vector 
connecting deuterons inside this molecule. Function V stands for the potential of the dfi-D2 
interaction |^ , leading to ddfi resonant formation. Vector Q2 connects the dfi and D2 centers 
of mass. We neglect contributions to the potential V from the molecules other than the 1-th 
molecule because we assume here that distances between different molecules in the target 
are much greater than the D2 size. The kinetic energy e of the dfi atom and its momentum k 
in the target frame are connected by the relation 

e = ky2Md^ . (11) 
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The Hamiltonian if of a pure D2 target, corresponding to the initial target energy £"0, 
has the form 



where Rj is the position of j-th molecule center of mass in the target frame (Fig. |^), U^y de- 
notes interaction between the j-th and j'-th molecule, and Mmoi is the mass of a single target 
molecule. 

The coordinate part \l/tot of the initial wave function of the system can be written as 
a product 

*tot = V'^^(ri) V^^f (^1) exp(2k ■ R,^) |0) , (13) 

where |0) stands for the initial wave function of the condensed D2 target, corresponding the 
total energy E^. Eigenf unctions of the operators B.^^ and E.Yi-2. are denoted by V'^^ and ip'^^', 
respectively. Using the relation R^^ = R; + the wave function \l/tot takes the form 

^tot = ^^^(ri) ^^f (ei) exp(2k ■ Q,) exp(^k ■ R,) |0) , (14) 

which is similar to that used in the case of dd^ formation on a single D2, except the factor 
exp(ik-R;) |0). This factor depends only on positions of mass centers of the target molecules. 

After formation of [{dd^)dee] complex, the total Hamiltonian of the system is well ap- 
proximated by the operator H^^^ 



'tot 



HU = H,,^{r, R) + (q) + V{q, r,R) + H , (15) 



where H^^^ is an internal Hamiltonian of dd/j, molecular ion, vectors r and R are its Jacobi 
coordinates. Relative motion of ddfi and d in the complex is described by a Hamiltonian 
which depends on the respective internal vector q. The final Hamiltonian H of the target, 
with the eigenfunction \ri) and energy eigenvalue En, is expressed by the formula 

2Mc ^'^^2M^oi " 

JT'l 3 3+3 (^i^ 

M^oA 1 _^2 +H = AH + H, 



Mc J 2M^oi • 

where Mc is the mass of the complex. The respective coordinate part "^[^^ of the total final 
wave function of the system is 

^;ot = ^i(r,R)C'^^(e)|^). (17) 

where ipj^fi i^^J^^ denote eigenfunctions of the Hamiltonians Hddii and H^, respectively. 

The energy-dependent resonant dd^ formation rate ^u!Ki,i'fKfi^) the condensed target, 
for the initial |0) and final \n) target states and a fixed dfi total spin F, is calculated using 
the formula 

^'S.,u,K,X^) = '^^NB^j \Aojn? - Sif + Eo- K) , (18) 
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with the resonance condition 

e + Eo = eif + En, (19) 

taking into account the initial and final energy of the target. The resonant energy for a free 
D2 is denoted by Eij and the transition matrix element is given by 

AiOJn ={Kt\V\^tot). (20) 

Using Eqs. (0) and ([TTD the matrix element (|20|) can be written as a product 

Aiojn = {n\ expitk ■ Ri)\0) %(£) (21) 

where Vif{e) is the transition matrix element calculated for a single D2 molecule The 
rate (|TBp can be additionally averaged over a distribution p„g of the initial target states at 
a given temperature T and summed over the final target states, which leads to 

>^u5<.,u,Kfie) = 2nNB,j\V,fis)\'J2pno |(n|exp(^k- R,|0)|2 

n,no (22) 
X 6{e - Eif + Eq- En). 

Factor Bif, defined by Eqs. (P), is due to the averaging over the initial projections and 
summation over the final projections of spin and rovibrational quantum numbers. This 
factor takes also into account a symmetrization of the total wave function of (i//+D2 system 
over three deuterium nuclei. 

Now we introduce a time variable t to eliminate the S function in the equation above and 
then we involve time-dependent operators, which is familiar in scattering theory (see, e.g., 
Refs p8| , p9|). Using the Fourier expansion of the 6 function 

6{£ -eif + Eo-En) = — dt exp(-it{e - e,f + Eo - E^)) (23) 

J -00 

one has 



/oo 
dt exp(-it {e - Sij)) p„o 

X (0| exp(— zk ■ Ri)|ra)(ra| exp{itEn) exp(ik ■ R;) exp(— zti?o)|0) • 

(24) 

Assuming that the perturbation operator AH is well-approximated by its mean value 

AH ^ {0\AH\0) =Ae,f = -{l~ M^jMc) < , (25) 

which is valid when the target relaxation time is much smaller than the dd^ lifetime of the 
order of 10~^ s, the matrix element in Eq. (^) can be expressed as 

(n| exp(zti?„) exp(zk ■ R^) exp(— it£'o)|0) 
= (n| exp {it{H + AH)) exp(ik ■ R;) exp(-zti7) |0) 
~ (n \ exp(ztA£:jj) exp{itH) exp(ik ■ R;) exp{—itH)\0) 
= {n\exp{itAeif) exp(2k- R;(t))|0) , 
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where R/(t) denotes the Heisenberg operator and £t in formula (^) is the mean kinetic 
energy of the target molecule at temperature T. 

Using the identity |^^)(^^| = 1 in Eq. (|2^ ) we obtain 



/oo 
dt ei,v(~it{e - e[f)) 
-oo 



(27) 

X (exp (-ik ■ (0)) exp (^k • R; (t)) )^ , 



where {■ ■ ■)t denotes both the quantum mechanical and the statistical averaging at temper- 
ature T, and e[j being the resonance energy 

e[f = e^f + ^e^f , (28) 

shifted by Ae^j < 0. Note that such a resonant energy shift was neglected in papers P7[ , 
where absorption of neutrons and 7-rays by heavy nuclei were considered. An estimation of 
the shift in the case of 7 emission from a nucleus bound in a solid, similar to Eq. (|25|) was 
given in Ref. |]30 |. 



A self pair correlation function Gs{r,t) is defined by the following equation 114 



(exp(-zk-R;(0))exp(zk- R(t)))^ = ^ ^(r, t) exp(ik ■ r) , (29) 

and the single-particle response function Si{K,u!) is given by the formula 

Si{K,u) = — I d^r dt Gs{r,t) exp(i{K ■ r — out)) . (30) 
2n J 



Thus, by virtue of Eqs. (^) and (pOf), the resonant formation rate in a condensed target can 
by expressed in terms of the response function 

= 27rNB,f\V.j{e)\'S,iK,u;) , (31) 

where the momentum transfer k, and energy transfer uj to the target are defined as follows 

K = k , u = e — e[j: . (32) 

The advantage of the Van Hove method is that all properties of the target, for given momen- 
tum and energy transfers, are contained in the factor iSj(k, cu). It is possible to rigorously 
calculate Si in the case of a perfect gas and in the case of a harmonic solid. However, a liquid 
target or a dense gas target is a difficult problem to solve. 

Proceeding as above one can obtain a similar formula for ^ufxf UiKi i^ ^ condensed target 
(in general, d^ spin F' after back decay can be different from d^ spin F before the formation) 

/d'^k ~ 
^|V^,(.)P5.(..,^'), (33) 

ijj' = e\j — e , e'ij = Eif + Aiif , 
Si is the response function calculated for the state \n) and 

Asif = {n\AH\n) = - (Mc/M^oi - 1) , (34) 
where St denotes the mean kinetic energy of the complex bound in the target. 
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IV. RESONANT FORMATION IN A HARMONIC SOLID 



It has been shown by Van Hove |TJ] that the self correlation function in the case of a gas 
or a solid with cubic symmetry takes the general form 



G,(r,t) 



3/2 



exp 



V 27(t) J 



(35) 



For a cubic Bravais lattice, in which each atom is at a center of inversion symmetry, 7(t) is 
given by the formula 



7(t) 



aw n [w) exp[—twt) , 

w 



(36) 



where Z{w) is the normalized vibrational density of states such that 

dw Z{w) = 1 , Z{w) = for W > t^max , 

Z(—w) = Z(w) , 



n„ (w) is the Bose factor 



n^{w) = [exp{f3w) - l]'' , f3 = (A;bT)-i . 



(37) 



(38) 



and the Boltzmann constant is denoted by /cb- 

The response function (pOD, after substitution of Eqs. (pSf), (|36|) and integration over r, 
can be written as follows 

(39) 



5i(K,^) exp ("-—^7(00)) 



X 



dt exp{—iut) exp 



( 



2M, 



mol 



[7(00) - 7(t)] 



7(00) denotes the limit of 7(t) at t ^ 00. This formula can be expanded in a power series 
of the momentum transfer n, which leads to 



Si{K,uj) = exp{-2W) 



{2wy 



n=l 



where 2W is the Debye- Waller factor, familiar in the theory of neutron scattering. 



2W 



K 



and the functions Qn are given by 

giiw,T) 



7(00 j 



2M 



mol Jo 



dw 



Ziw) 



w 



coth {\(3w) 



(40) 



(41) 



1 Ziw) 



KH + 1] 



7(00) w 

POO 

gn{w,T)= / dw' gi{w - w\T) g.a-i{w' ,T) 



(42) 



dwgn{w) = 1 . 



In the case of a cubic crystal structure 2W can also be expressed as 

2W = i{0|u^|0)K^ 



where u is the displacement of a molecule from its lattice site. Substitution of Eq. 
Eq. ( piD leads to the following formation rate 



2TxNB,f\V,f{e)\^e^^{-2W) 



1 



{2W)' 



n- 



(43) 

to 

(44) 



The first term in expansion (|4^) represents a sharp peak describing the 6 profile recoilless 
formation. The next terms give broad distributions corresponding to subsequent multi- 
phonon processes. In particular, the term with n = 1 describes formation connected with 
creation or annihilation of one phonon. 



If 2W ^ 1 we deal with so-called strong binding p6| where only the few lowest terms in 
the above expansion are important. On the other hand, in the limit 2W ^ 1 (weak-binding) 
many multi-phonon terms give comparable contributions to (pj). Therefore, for sufficiently 
large it is convenient to use the impulse approximation in which j{t) is replaced by its 
value near t = 



7(t) ~ 7(0) + tt 3 
This leads to the asymptotic formula for Si 

1 



2 p 



where 



A 



exp 



n 



2M, 



mol 



(45) 



(46) 



(47) 



The mean kinetic energy £t of a molecule in the solid, which also determines the resonance 
energy shift (p5|), is equal to 



dw Z{w) w {w) + |] . 



(48) 



The energy £t contains a contribution from the zero-point vibrations and it approaches 
"ik-oT /2 only at high temperatures T ^ Wmax/^B- Function (|4^) is a Gaussian with re- 
sponse centered at the recoil energy TZ. Therefore in the weak binding region the resonant 
formation rate takes the Doppler form obtained by Bethe and Placzek ^ for resonant ab- 
sorption of neutrons in gas targets ||31|. However, the resonance width (^7]) in the solid at 
temperature T is different from the Doppler width in a Maxwellian gas Agas = 2y/kBTTZ 
unless the temperature is sufficiently high. This phenomenon was pointed out by Lamb in 
his paper concerning resonant neutron absorption in solid crystals. By virtue of the 
equations above one can introduce for the solid an effective temperature Tgfj 



i-eS 



B • 



(49) 



^ In fact, formula 
width r 0. 



3) is the hmit of the Bethe formula in the case of a very narrow natural resonance 
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V. RESONANT FORMATION IN FROZEN DEUTERIUM 



The following considerations concern the solid deuterium crystals used in the TRIUMF 
experiments [0, Q , though the results presented below can be applied to targets obtained 



in similar conditions |12, 34 1. At TRIUMF thin solid deuterium layers have been formed by 



rapid freezing of gaseous D2 on gold foils at T = 3 K and zero pressure. According to Ref. |3^ 
such deuterium layers have the face-centered cubic (fee) polycrystalline structure. Since the 
distance between the neighboring molecules is a few times greater than the diameter of 
a D2 molecule and the Van der Waals force that binds the solid is weak, one can neglect 
perturbations of the resonant formation potential V due to these neighbors. 

The deuterium crystals at zero pressure are quantum molecular crystals. The amplitude 
of zero-point vibration at 3 K equals 15% of the nearest neighbor distance. A single-particle 
potential in this case is not harmonic and the standard lattice dynamics leads to imaginary 
phonon frequencies. However, the standard dynamics can be applied after a renormalization 
of the interaction potential, taking into account the short-range pair correlations between 
movement of the neighbors ||3^. In result, the theoretical calculations |^ of the phonon 



dispersion relations give a good agreement with the neutron scattering experiments |^ and 
the Debye model for solid deuterium can be used as a good approximation of the phonon 
energy distribution 

10 if w > , 

with the Debye energy = /cbQd and Debye temperature Gd taken from the neutron 
experiments. For T = 3 K we use the Debye model of an isotropic solid with 6d = 108 K 
corresponding to the maximal phonon energy = 9.3 meV. Thus, we are dealing with the 
limit T/Od <^ 1 where 

7(oo) = |u;-\ £t = ^w^^5.2meV, Teg = | 60 ~ 40 K , (51) 



are very good approximations of Eqs. (^l]), ( |48| ) and (|49D. The Debye- Waller factor and 
mean kinetic energy St at lowest temperatures are determined by contributions from the 
zero-point D2 vibration in the lattice, and therefore these quantities do not tend to zero at 
T — > 0. The zero-point energy is not accessible energy, but its effects are always present. 

The values of the resonance energies depend on initial and final rovibrational quantum 
numbers of the system. In solid hydrogens at low pressures these quantum numbers remain 
good quantum numbers, but excited energy levels broaden into energy bands (rotons and 
vibrons) due to coupling between neighboring molecules |TB[. The calculations presented in 
the literature concern pure solid H2, HD and D2 targets and only lowest quantum numbers. 
The problem of a heavier impurity, such as {dd^)d complex in D2, has not been considered 



yet. However, knowing that the width of the rotational bands can reach about 1 meV [|18 
a possible influence of this effect on the calculated formation rates and fusion neutron time 
spectra is discussed in the next section. 

At low temperatures all D2 molecules are in the ground vibrational state z/j = and dd/j; 
is formed via the excitation of the complex to the state Uf = 7. Unless a catalyst is applied, 
rapidly frozen deuterium is a mixture of ortho-D2 {Ki = 0) and para-D2 {Ki = 1). In the 
TRIUMF experiments gaseous deuterium was pumped through a hot palladium filter before 
freezing. Therefore the solid target was a statistical mixture (2:1) of the ortho- and para- 
states {Ki = stat). Since the para-ortho relaxation without a catalyst is very slow (0.06%/h) 
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in solid deuterium the population of these states is not changed during experiments of 
a few days. 

The lowest resonance energies Eif and e[j, for fixed Ui, Uf and different values of F, Ki, 
S and Kf are shown in Table | |jl0[. A few of them have negative values, which means that 
to satisfy the resonance condition e = Eif an energy excess in the dfi+D2 system should be 
transferred to external degrees of freedom. This is possible in dense targets, where energy 
of neighboring molecules can be increased. Such an effect, due to triple collisions in gas 
targets, has been firstly discussed in Ref. In a solid, the energy excess is lost through 
incoherent phonon creation. According to (P3|), (^), and (^Tj), in the considered 3 K solid 
deuterium all resonant energies e[j are shifted by Asif ~ —1.81 meV. One can see that 
all resonances for F = ^ are placed at higher energies, which is caused by dfi hyperfine 
splitting AE^^^ = 48.5 meV. All resonance energies £-j < ~ 10 meV are connected with 
formation from the upper spin state F = | of dfi. However, only resonances corresponding 
to the dipole transitions Ki = Kf = 1 and i^'j = 1 — > Kf = 0,2 can give a significant 
contribution to the formation rate at lowest energies. Other transition matrix elements 
described in Ref. tend to zero when e —>■ (see Figs. ^ and ^ obtained for Ki = and 



K, = l). 

The low energy rates {e < wd) are calculated using formula (0) with a few most sig- 
nificant terms of the response function expansion (^) taken into account. Fig. ^ shows the 
function Si{K,e — e'if) corresponding to the two dipole transitions in para-D2. The sub- 
threshold resonance, with e[f k, —9.0 meV, gives contributions to the formation rate only 
through the phonon creation processes. For e^f k, l.Q meV, the non-phonon process is pos- 
sible and it is represented by a vertical line. Different peaks in this figure describe processes 
connected with different numbers of created phonons. In particular, one-phonon processes, 
which are proportional to Z{w) with the characteristic Debye cutoff, can be clearly distin- 
guished. Since the n-phonon term in ( ^OD is proportional to k^", the dd^ formation rate 
tends to zero at e — 0. Note that the phonon annihilation gives negligible contribution to 
the rate at very low target temperatures T <^ 6d. 

In order to compare the calculated formation rates with experiments the summed rates 

^KM) = Y.^'ku,Kr ^. = 0, z.; = 7. (52) 



XkM) introduced 



3 

2- 



In Fig. ^ the formation rates A^^ {e) in the solid ortho-D2 and para-D2 are shown for F 
In the case of resonances satisfying the condition e^f < we have 2W < 1 and the 
expansion ( ^^ is used. The two strong peaks represent the recoilless formation process, 
without phonon excitations. The delta function profile of every peak is shown as a rectangle 
with a height equal to the formation rate strength divided by the total decay width (^ 
0.8 X 10~^ meV). The strength defined as the value of the factor standing before S{uj) in 
the expansion (PD, is equal to 0.1061 eV-/is^^ for the resonance Ki = ^ Kf = 1 in 
solid ortho-D2. The transition Ki = 1 —* Kf = 2 in para-D2 gives 0.07544 eV-/is~^ as 
the resonance strength. Higher resonance energies involve many multi-phonon terms and 
therefore we use the asymptotic form (|46|) of Si for e^f > . All formation rates presented 
in the figures are normalized to the liquid hydrogen density Nq = 4.25 x 10^^ atoms/cm^. 

Though in Monte Carlo simulations, involving energy-dependent rates of different pro- 
cesses, the "absolute" formation rates (e) should be used, it is convenient to introduce an 
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effective formation rate A^. (e) which leads to the nuclear dd fusion in [{dd^)dee\ complex. 
Back decay of the complex to the (i/i+D2 system, characterized by the quantum numbers K[ 
and F', strongly influences the fusion process because the back-decay rates are comparable 
with the effective fusion rate A/ ~ 374 /is~^ 0. Since in a solid target rotational deexci- 
tation of the asymmetric complex is much faster than back decay and fusion, it is assumed 
that back decay takes place only from the state Kf = 0. The effective formation rate is then 
defined by the following formula 



Kf,S 



fus 



where the fusion fraction is given by 



i/i = 0, i^f = 7, 



(53) 



fus 



A/ 



= Ay + J] r^^' 



E 



i UfKf,u,Ki 



(54) 



Since the frequency of lattice vibrations (~ w^/h ~ 10'' /is~^) is many orders of mag- 
nitude greater than the back-decay and fusion rates, energetic phonons created during 
the ddfi formation process are dissipated. At 3 K the number of phonons with energies 



w > k-gT 0.26 meV is strongly suppressed by the Bose factor 



Therefore back 



decay with phonon annihilation at T <C is negligible. In particular, the phonon channel 
of decay of ddfi, formed from rf/i state = | due to the subthreshold resonances, is closed 
because this would require an annihilation of a phonon with energy of a few meV. In this case 
back decay is connected with the spin-flip transition to F' = i. Since the corresponding en- 
ergy release of a few tens of meV is much greater than the Debye energy [AE^^^ ^ w/)), the 
ddfj, decay rate is dominated by contributions from simultaneous phonon creation processes. 
After integration of formula (^3]) over direction of vector k one obtains 



ySF' 

VfKf,ViKi 



A 



if 



dk e\V^f{e)\''Si{k\u}'^ 



(55) 



and then substitution of expansion ( ^OD and integration of the recoilless term lead to 



ySF' ^ A/ 



+ Y, r dkk'\V,jis)\'expi-2W)gr.iu;',T) 



n=l 



(56) 



where 



2W 



2M, 



c 



2W, 



if 



2W{k, 



if) 



if 



2Me' 



if 



(57) 



It is assumed in the formula above that the phonon energy spectrum of solid deuterium 
containing [{ddn)dee] is similar to that of a pure deuterium lattice. The problem of lattice 
dynamics of a quantum solid deuterium crystal containing a small admixture of a heavier 
isotope has not been considered yet in literature, at least to the knowledge of the authors. 
However, this approximation is reasonable since the Debye temperatures of solid hydrogen 
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and deuterium at 3 K are very similar independently of the mass difference of these 



isotopes. Therefore it is assumed that during the ddfi lifetime the mean kinetic energy St of 
the complex reaches the energy St characterizing a pure deuterium solid. Thus the resonance 
energy shift (133) is approximated by 



Ae^f ^ - (Mc/M^oi - 1) ~ -2.77 meV , (58) 

which gives e[jr = £j/ — 2.77 meV. 

The effective formation rates in 3 K solid deuterium for -F = | are shown in Fig. |^. 
The phonon part of the rates below a few meV is about two orders of magnitude lower 
than the average rate of 2.7 /is~^ derived from the experiment [0, |T^. This means that 
at e <^ Wd the phonon contribution to the total resonant formation rate is even smaller 
than the non resonant dd/i formation rate Anr ~ 0.44 /is"^ and that the estimation of the 
phonon contribution given in Ref. is strongly overestimated. Therefore, the experimental 
results can only be explained by resonant ddfi formation at energies s > 1 meV, where 
the rate exceeds significantly the value of 1 /xs~^. A cusp at 0.3 meV in para-D2 is due 
to the formation with simultaneous one-phonon creation, connected with the subthreshold 
resonance Ki = 1 Kf = 0. This implies a significant difference between the resonant 
formation in ortho-D2 and para-D2 below 1 meV. However, this difference is difficult to 
measure because of a broad distribution of dfi energy. Note that a similar subthreshold 
phonon effect in the case of resonant dt/i formation in solid deuterium has been discussed 



in Ref 41 



In the solid target the fusion fraction V^^^ ~ 0.3 and the total resonance width a; 
0.8 X 10-3 meV for both 5 = i and ^ = f . The back-decay rate T^^' from S = ^ to F' = ^ 
equals about 843 /is"^. Decay S" = | — F' = | is impossible. In the case of = | we have 
obtained T^^' ^ 281 /is"^ for F' = ^ and T^^' «i 610 /is"^ for F' = |. Phonon creation 
processes give dominant contributions to the back-decay rates, e.g., the non-phonon part of 
, given by the first term of expansion (^), equals 169 fis~^. Therefore the dfi energy 
spectrum, after back decay in the solid, is not discrete. 

In Fig. H the effective rates in solid deuterium for F = ^ are presented. For the sake 
of comparison the formation rate for 3 K ortho-D2 gas is also plotted. The "gas" curve 
has been calculated using the asymptotic formula ( ^61) for Si with T^s = 3 K.This figure 
shows that in a real solid deuterium target the rates are smeared much more than in a gas 
target with the same temperature, because of the zero-point vibrations. Therefore even at 
relatively high dfi energies of some 0.1 eV one should not neglect the solid effects and use 
the formation rates calculated for a 3 K Maxwellian gas. 



VI. MONTE CARLO CALCULATIONS 

The calculated energy- dependent ddfi formation rates have been applied in our Monte 
Carlo simulations of fiCF in 3 K solid deuterium targets. The final dfi energy distribution 
after back decay, including simultaneous phonon creation processes, has been determined 
through a numerical integration of Eq. (^61). The calculated distribution is shown in Fig. |^ 
for S" = |, Kf = and F' = |. The rotational transitions to K'^ = 0, 1, 2 with no phonon 
creation are seen as the delta peaks. The continuous energy spectrum describes phonon 
creation contribution to dfi energy. Note that, opposite to ddfi formation rates, this phonon 
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contribution (for a given rotational transition peak) extends towards lower energies. The 
average dfi energy after dd^ back decay equals about 30 meV, for the presented spectrum. 

The dd fusion neutron and proton spectra depend on the time evolution of dfi energy. 
This energy is determined by differential cross sections of different scattering processes of dfi 
atoms in a given solid target, including elastic scattering, rovibrational transitions, spin-flip 
reactions and phonon processes. The scattering cross sections in a solid are calculated using 
the Van Hove method. Some results of such calculations for dfi atoms in fee solid deuterium 
have been presented and discussed in Ref . . The incoherent processes, such as spin- flip or 
rovibrational transitions, are described by the self pair correlation function Gs{r,t) defined 
by Eq. (^) . The Bragg scattering and coherent phonon scattering are connected with a pair 
correlation function G{r,t) |]14| . 

In Fig. |10| is shown the total cross section for dfi{F = |) scattering in the statistical 
mixture of 3 K solid ortho-D2 and para-D2. Bragg scattering, with the Bragg cutoff at 
£3 = 1.1 meV, and incoherent elastic scattering do not change d/i energy because of the very 
large mass of the considered solid target. Below 1.7 meV the d^ atom is effectively acceler- 
ated, mainly due to the rotational deexcitation of para-D2 molecules 0. This transition 
is enabled by muon exchange between deuterons in diJ,+D2 scattering. The curve "0 1" 
in Fig. |10|, describing the rotational deexcitation, includes contributions from simultaneous 
incoherent phonon processes. This cross section at e = 2.5 meV equals 0.22x10"^'' cm^, 
which is about three times less (taking into account the statistical factor of 1/3 for K = 1 



states) than the estimation given in paper I^Tj. Phonon annihilation is a much weaker dfj, 



acceleration mechanism than the rotational deexcitation. 

Since the coherent amplitude for dfi elastic scattering on a single D2 molecule is greater 
by two orders of magnitude than the incoherent amplitude, the coherent processes involving 
conservation of momentum dominate low energy dfi scattering in solid deuterium. It is 
especially important below a few meV, where the shapes of coherent and incoherent cross 
sections differ strongly. The small phonon creation cross section below 1.1 meV, leading to 
dfj, energy loss, is due to the incoherent amplitude. Coherent phonon creation is impossible 
below e^. This limit is obtained in the case of coherent one-phonon creation process, for 
the total momentum conservation involving the smallest (non-zero) inverse lattice vector r, 
which also fixes the position of the first peak of the Bragg scattering at = 1.1 meV. For 
T = one-phonon creation is possible only if the dfi velocity is not lower than the sound 
velocity in the crystal, which is well-known in neutron physics. According to Ref. the 
mean sound velocity in solid deuterium equals about 1.2x10^ cm/s and this corresponds 
to dfi energy of 15 meV. Therefore, neglecting the inverse lattice contribution to the one- 
phonon creation cross section in Ref. ^l|] leads to the severe underestimation of dfi slowing 
down at lowest energies and subsequent overestimation of dfi kinetic energy. 

Above 1.7 meV phonon creation already prevails over all acceleration processes. However, 
the effective deceleration rate below is strongly suppressed by the dominating Bragg 
elastic scattering. At energies above some 10 meV subsequent rotational and then vibrational 
excitations of D2 molecules become important and they provide a very fast mechanism of 
dfi deceleration at higher energies. 

The total cross section for dfi{F = |) scattering in a pure 3 K ortho-D2 target (see 



Fig. |TT]) is quite similar to that shown in Fig. |T0[ A significant difference is the lack of 
rotational deexcitation. Therefore phonon annihilation is the only, and weak, acceleration 
mechanism. It dominates the inelastic cross section below 1.4 meV. 

Fig. |T2| presents the time evolution of average dfi{F = |) atom energy ^avg, obtained from 
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our Monte Carlo calculations. It is assumed that the target is infinite and that djj, atoms 
have initially a Maxwellian energy distribution with a mean energy of 1 eV. A statistical 
initial population of dfi total spin is used and the theoretical non-resonant part of the total 
spin-fiip rate As i is multiplied by a single scaling factor of 0.4, in order to keep agreement 
with the experimental values |10|, ^ of the spin-flip rate. The calculations have been 
performed for ortho-D2, para-D2 and their statistical mixture (stat). One can see that 
dn mean energy of 10 meV is reached already after 5 ns. Then, below the Debye energy, 
deceleration become very slow. The lowest value of ^avg is determined by the intersection 
point of the cross sections of the acceleration processes and phonon creation process. In the 
case of a statistical mixture ^avg ~ 1-7 meV, for = we have ^avg ~ 1-4 meV. Finally, for 
pure para-D2, with a contribution to the total cross section from the rotational transition 
K = 1 ^ three times greater than that shown in Fig. |10|, ^avg ~ 2.2 meV. Thus, d^ atoms 
are never thermalized and their energy is signiflcantly greater than 1 meV. For para-D2 the 
mean energy is always greater than the energy of the lowest resonance peak e[j = 1.6 meV. 
However, even if ^avg is smaller than e[j-, a signiflcant part of dfi atoms has energy e > e[f 
because of a large admixture of hot d^ atoms at t = |l^, |16| and slow deceleration below 
10 meV. 

Since at energies of a few meV the lowest delta peaks are dominant in the resonant 
formation, their contributions to the mean effective formation rate are shown in Fig. |l^ for 
gas and solid deuterium (stat) targets, assuming steady Maxwell distributions of d^{F = |) 
energy, with different ^avg- The maximum average rate of about 6 fis~^ in the solid is 
due to the resonance energy shift of —1.8 meV. The experimental result of 3 /is~^ can be 
explained because ^avg is greater than 1 meV. However, in order to obtain large fusion 
neutron and proton yields through resonant ddfi formation, the width F'^ of the resonance 
peaks in solid can not be too narrow. The peak resonant rates of a few 10^/is~^ have been 
obtained assuming the discrete values of the rovibrational D2 energies in solid deuterium 
and F*^ ~ 10~^ meV. These resonant rates are many orders of magnitude greater than the 
inelastic scattering rate ~ 10 /is~^. In such a case dfi atoms are very quickly (compared to 
dfi{F = |) lifetime) removed from the regions of resonance peaks and the contribution of the 
recoilless resonances to the neutron yield is negligible. The Monte Carlo simulations have 
shown that the neutron yield from the phonon part of the resonant rates gives only some 10% 
of the yield observed in the experiments. In result, the calculated time spectra, obtained 
for the small F"^, are dominated by weak non-resonant ddfi formation, which disagrees with 
the experimental data. Therefore, we have investigated influence of a broadening of the 
non-phonon resonant peaks, due to the presence of molecular rovibrational bands in solid, 
discussed in Ref. [|18[. Since in the literature there is no information concerning the proflle 
of such bands, we have assumed a rectangular shape of the resonance peaks. The resonance 
strengths have been flxed and their widths have been changed in the limits 0.001-1 meV. It 
turns out that good Monte Carlo results are obtained for F'^ ^ 0.5 meV, which is consistent 
with the rotational bandwidths of about 1 meV reported in Ref. ||18[. This gives the resonant 



formation rate of 294 fis~ for the recoilless peak in ortho-D2, and respectively 214 /is~ in 
para-D2. In Fig. |1^ one sees the resonant formation rate at lowest energies for F"^ = 0.5 meV 
and for the statistical mixture of ortho- and para-states. Also shown is the Monte Carlo 
distribution of dfj,{F = |) energy, calculated for times t = 10 ns and t = 30 ns. The 
Maxwell distribution of initial dn energy, with ^avg = 1 eV, has been assumed. Two minima 
in the dn energy distribution appear quickly at the positions of the resonance peaks since 
the respective ddfi formation rates are comparable with the total inelastic scattering rate of 
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about 30 /is~^. 

The dd fusion neutron spectrum, calculated assuming the same initial d^ energy and 
resonance profiles, is shown in Fig. A 3.2x10"^ concentration of nitrogen is included in 
order to fit the TRIUMF target conditions. The solid line plotted in this figure has been cal- 
culated using the steady-state kinetics model with the effective formation rate A^jfj = 3 /is^^ 
and total spin-flip rate As i = 36 /is~^ taken from the fits to the experimental data ||13[. The 



2 ' 2 



slope of the spectrum at t < 80 ns is determined by the rates XJtat^ "^§.5' '^^ scattering 
rate which also changes the population of dn^F = |) atoms in the vicinity of the resonant 
peaks. The steady-state kinetics model does not include the process of dfi deceleration. 
Therefore, fits using this model could entangle the deceleration rate with the formation and 
spin-flip rates. The mean formation rate, calculated directly in the Monte Carlo runs, is 
a function of time, and it stays at the level of 1-3 /is~^. The spectrum slope at large times 
t > 100 ns, when dfi{F = |) atoms practically disappear, are due to the nonresonant dd/i 
formation from = | and to the muon transfer to nitrogen contamination. 

The shape of the time spectra practically does not change when the mean energy ^avg of 
the initial single Maxwell distribution varies in the limits 0.01-1 eV. On the other hand, the 
spectra change strongly if a significant part of dfi atoms at t = has energy smaller than 
the energy of the lowest resonant peak, which can be observed using a more complicated 
(e.g. two-Maxwell distribution). Assuming that F'^ is greater than 0.5 meV we obtain results 
which begin to differ significantly from the analytical curve calculated with the experimental 
parameters. In particular, the ratios of neutron yields from the short and large times begin 
to disagree. Fits of the calculated spectra to the experimental data would enable a better 
determining of F'^ and a shape of the initial dfi energy. However, this is not the purpose 
of this work. A qualitative comparison of Monte Carlo spectra with the experimental data 



has already been performed in article ||T3[. In this case good fits were not obtained since 



at that time the resonant ddfi formation rates in solid D2 and dfi scattering rates including 
coherent effects in the solid were not yet available. 

Our calculations show that strong resonant ddfi formation takes place both in ortho-D2 
and para-D2. There are certain differences between the neutron time spectra from these 
targets (see Fig. |l^), caused by the different positions and strengths of the lowest resonance 
peaks. Also dfi slowing down process differs slightly in the two cases. The neutron yield at 
larger times is smaller for ortho-D2 since in this case the resonance peak is placed at higher 
energy of 2.3 meV. Therefore, dfi atoms are removed faster from the peak compared to the 
situation in para-D2, where the resonance is observed at 1.6 meV. A greater mean dfi energy 
in para-D2 (cf. Fig. 0) leads also to a stronger overlap of the resonance peak and dfi energy 
distribution at t > 20 ns. However, the differences between the spectra can be clearly seen 
only in high-statistics experiments. 



VII. CONCLUSIONS 

The methods used for description of resonant neutron and 7-ray absorption in condensed 
matter have been directly applied for calculation of resonant ddfi formation and back-decay 
rates in condensed deuterium targets. These rates are expressed in terms of the Van Hove 
single-particle function, which depends on properties of a given target. In particular, we 
have derived the analytical formulas for the rate in the case of resonant ddfi formation in 
a harmonic solid deuterium. The calculations show great differences between resonant ddfi 
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formation in 3 K solid deuterium and in 3 K D2 gas. In solid, the formation at a few meV, 
which determines the experimental results, is dominated by presence of the strong recoilless 
resonant peaks. On the other hand, the formation with simultaneous phonon creation is 
important above the Debye energy. The resonance profiles in the solid at higher energies are 
similar to that in D2 gas, but with the effective temperature equal to 40 K. This temperature 
is determined by the energy of zero-point vibration of D2 molecule in the lattice. Phonon 
creation is always important in the case of ddfi back decay because it is connected with 
energy release of a few tens meV, which is much greater than the Debye energy. 

A condition T/G^ ^ 1 is fulfilled for any solid deuterium target at low pressure. There- 
fore, the parameters determining solid state effects (Debye- Waller factor, mean energy of D2 
vibration in solid) weakly depend on target temperature T. They are expressed in terms of 
the Debye energy w^^ which does not significantly change with the varying solid tempera- 
ture T. In result, the resonant ddfi formation rates in solid deuterium for different T are very 
similar and one may expect that the average formation rates, derived from measurements 
performed at different temperatures, will also be very close. This is confirmed by the results 
of experiments carried out at TRIUMF and at JINR. 

The structure of a solid deuterium target depends on its temperature and history. Targets 
maintained at T > 4 K have the hep structure [^. Though our calculations have been 



performed for fee crystals, the obtained results are also good approximations of the resonant 
rates in hep polycrystals since the Debye temperature and nearest neighbor distance are 
similar for these two lattices. In general, the formulas derived in this paper can be used in 
a wide range of target temperature and density, with appropriate experimental values of the 
Debye temperature and lattice constant taken into account. 

The Monte Carlo calculations show that dfi deceleration below the Debye energy is very 
slow and that mean energy of dfi{F = |) atom is always significantly greater than 1 meV. The 
energy distribution of (i/i's during their lifetime is very broad (at least a few meV), therefore 
a strong overlap of this distribution and lowest resonance peaks takes place, leading to a large 
mean ddn formation rate in solid deuterium. However, explanation of the experiments is 
possible only if the broadening of rovibrational molecular levels in solid is taken into account. 
We obtained reasonable results assuming that the strengths of the recoilless resonant peaks 
are constant and that the rotational bands increase the resonance peak width to 0.5 meV. 



Note that, according to Ref. |T8[, high pressures lead to a greater broadening and even to 
a mixing of rotational states. This could complicate a comparison of theory and high-pressure 
experiments. The phonon part of the resonant rate give only about 10% contribution to the 
calculated neutron time spectra. 

The dd fusion neutron spectra calculated for ortho-D2 and para-D2 solid targets are 
quite similar. Small differences between the spectra are due to the different energies and 
strengths of the lowest resonant peaks, and to a slightly higher mean d/i energy in para-D2. 
These differences can be clearly seen only in high-statistics experiments. Our calculations 
do not confirm a lack of strong resonant ddn formation in solid ortho-D2, predicted in the 



papers In order to verify the theory it is necessary to perform measurements in 

pure ortho-D2 and para-D2 solid targets under the same conditions. 
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FIG. 1: System of coordinates used for the calculation of resonant formation of the com- 
plex \(ddii)dee] in a condensed deuterium target. 




FIG. 2: Position of impinging dji atom with respect to the condensed target. 
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FIG. 3: Transition matrix elements |V^/(£)p for Ki = and Kf = 0, 1, 2 versus d/x energy 
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FIG. 4: Transition matrix elements for = 1 and Kf = 0, 1, 2 versus d/x energy 
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FIG. 5: Response function Si{K,e — e[j) (in arbitrary units) for the para-D2 crystal at 3 K. The 
dashed hne is obtained for the subthreshold resonance e'^ ~ —9.0 meV, the soHd hne corresponds 
to e'^f K. 1.6 meV. The vertical line represents the rigid lattice term 5{e — e'-f) exp(— 




FIG. 6: Formation rate A^. (e) for F = | in 3 K ortho-D2 (solid line) and para-D2 (dashed line). 
The labels "1 2" and "0 — > 1" denote the rotational transition Ki — > Kf corresponding to the 
lowest non-phonon processes. 
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FIG. 7: Effective formation rate (e) for F = | in 3 K solid ortho-D2 and para-D2. The labels 
"1 — > 2" and "0 — > 1" denote the rotational transition Ki Kf corresponding to the lowest 
non-phonon processes. 




FIG. 8: Effective formation rate A^. (e) for F = ^ in 3 K solid ortho-D2 and para-D2. The 
label "gas" denotes the curve obtained for 3 K gaseous deuterium {Ki = 0), using the asymptotic 
formula (|4^ ) for the response function Si with T^g = 3 K. 
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FIG. 9: Distribution of final dfi energy after ddfi back decay from S = ^,Kf = to F' = ^, 
K'- = 0, 1,2. The three peaks describe the rotational transitions without a simultaneous phonon 
excitation. 
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FIG. 10: Total cross section for d/x(F = |) scattering in statistical mixture of solid ortho-D2 
and para-D2. The label "1 0" denotes the rotational deexcitation = — > 1 of a target D2 
molecule. The curves "—phonon" and "+phonon" stand for scattering with phonon annihilation 
and creation, respectively. The Bragg cross section is calculated for the fee polycrystalline lattice. 
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FIG. 11: Cross section for dfi{F = |) scattering in solid ortho-D2. The labels are identical to those 



in Fig. 10 
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FIG. 12: Calculated time evolution of average dfj, energy e^vg for F = | in 3 K solid ortho-D2, 
para-D2, and their statistical mixture (stat). A Maxwell distribution of dfj, initial energy, with 
mean energy of 1 eV, has been assumed. 
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FIG. 13: The effective resonant ddfj, formation rate as a function of mean CMS energy ^avg of 
d^i{F = |) atom for gas and solid deuterium targets. A steady Maxwell distribution of d/i energy 
is assumed for a given £avg- The contributions from the two lowest resonant peaks to the formation 
rate are taken into account. 
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FIG. 14: Resonant ddfi formation rate for F = | in the statistical mixture of ortho-D2 and para-D2 
for the resonance peak width F"^ = 0.5 meV. Monte Carlo distribution of dfx energy at t=10 ns 
and i=30 ns after the muon stop is plotted (in arbitrary units). 
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FIG. 15: The Monte Carlo fusion neutron spectrum for the statistical mixture of 3 K solid ortho-D2 
and para-D2 (solid line). The dashed line represents the spectrum obtained using an analytical 

~3/2 _i 

steady state kinetics model with A^^^^ = 3 /US . The initial dfi energy is given by a Maxwell 
distribution with mean energy of 1 eV. The width T"^ of the non-phonon resonances is fixed at 
0.5 meV. A 3.2xl0~^ concentration of nitrogen is included. 




FIG. 16: Calculated neutron spectra from 3 K solid ortho-D2 and para-D2. The Maxwell distri- 
bution of initial dfi energy with £avg = 1 eV and F'^ =0.5 meV have been assumed for the both 
targets. 
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TABLE I: The lowest resonance energies of ddfx formation in dfi scattering from single D2 
molecule (ej/) and from 3 K solid deuterium target (e^j)- These energies are given in the re- 
spective CMS systems. 

Eif (meV) (meV) F Kf S 
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